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Abstract 

A correspondence between the sextic anharmonic oscillator and a pair of 
third-order ordinary differential equations is used to investigate the phe- 
nomenon of quasi-exact solvability for eigenvalue problems involving dif- 
ferential operators with order greater than two. In particular, links with 
Bender-Dunne polynomials and resonances between independent solutions 
are observed for certain second-order cases, and extended to the higher-order 
problems. 
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1 Introduction 

Despite the simplicity of one-dimensional quantum-mechanical systems, full solv- 
ability is very much an exception rather than the rule. The archetypal example is 
the harmonic oscillator, for which exact solvability breaks down completely under 
almost any kind of perturbation. Nevertheless, Turbiner [T] found the surpris- 
ing fact that in some cases, corresponding to certain multi-parameter families of 
second-order differential eigenvalue problems, there are regions of the parameter 
space for which a finite subset of the spectrum can be found algebraically. Tur- 
biner and Ushveridze [2] dubbed these models quasi-exactly solvable, or QES. For 
the QES problems of [I], and almost all other examples, the differential operators 
% act invariantly in a finite dimensional subspace V n +\ spanned, possibly after a 
gauge transformation, by the monomials 1, x, . . . , x 11 . 

It is natural to ask whether similar ideas might apply to eigenvalue prob- 
lems involving higher-order differential operators. Turbiner has given a general 
classification of differential operators of order k with finite dimensional subspace 
'Pn+i [3] for which n + 1 eigenfunctions have the form of a polynomial of order at 
most n. Examples of higher-order QES operators arising in multi-mode bosonic 
Hamiltonians relevant to nonlinear optics may be found in [H [5j EJ El [8] . 

In this paper we consider third-order problems, motivated by a link between a 
specific family of third-order differential equations and the Schrodinger equation 
for the sextic anharmonic oscillator that was discovered and explored in [91 [12] . 
By tuning the parameters so that the second-order eigenproblem is at the points 
at which its spectrum is QES, the corresponding isospectral third-order eigen- 
problem must also be QES. It turns out that the quasi-exact solvability of these 
third-order differential equations is rather subtle, and even though a finite subset 
of eigenvalues can be found exactly, the most natural ansatz for the corresponding 
eigenfunctions turns out not to hold. Nevertheless, we shall exhibit a number of 
interesting properties of these and related second-order eigenproblems, and show 
in particular that the QES eigenfunctions are the generating functions for sets of 
Bender-Dunne polynomials whose zeros correspond to the QES eigenvalues. 

2 Isospectral second and third order eigenproblems 

Turbiner showed [ij that the sextic potential 



with regular boundary conditions imposed on the positive real axis at x = and 
x — > oo by requiring for I > —1/2 



^Ix-M) = x l+ \l + 0{x 2 )) , iP(x) = 0{x- z / 2 ~ a ' 2 e~ xi/A ) as x -> oo , (2.2) 





is quasi-exactly solvable along the lines 



a = aj = -{21 + 1 + 4J) , J = 1, 2, 3, . . . . 
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The operator 7i2(aj, I) is quasi-exactly solvable because it acts invariantly in the 
finite dimensional subspace (/o, /i, . . . , fj-x) where f n (x) = x l+1 exp(— x 4 /4) x 2n . 
The gauge transformation x _i_1 exp(x 4 /4) %2 x l+l exp(— x 4 /4) and variable change 
x 2 = w transforms the Schrodinger equation (|2.ip into a second-order differential 
equation which acts invariantly on Vj = (l,w, pQ. 

A convenient way to handle the quasi-exact solvability of (|2.ip is through the 
Bender-Dunne polynomials introduced in [10]. The idea is to write a solution to 
(|2.ip in the following factorised form 

n=0 

For f|2.4j) to solve (|2.ip . the coefficients P n must satisfy the recursion relation 

P n (E)=EP n „ 1 (E) + 16(n-l)(n-j-l)(n + l-l/2)P n „ 2 (E), (n > 1) (2.5) 

with j = j(a,Z) = -(a+2Z+l)/4 and P (£) = 1- Fr om ([23]), Pi = E, and P n is 
a polynomial of degree n in E 1 , known as a Bender-Dunne polynomial. As long 
as I — n— 3/2 for any n E Z + , f|2.4[) will yield an everywhere-convergent series 
solution to (|2.1|) . This solution automatically satisfies the boundary condition 
at the origin, but at general values of E, it will grow exponentially as x — > oo. 
However, if a and I are such that j(a, I) = J is a positive integer, the second term 
on the RHS of (|2,5p vanishes at n = J+l, and all subsequent P n factorise: 



P n+ j(E, a j, I) = Pj(E, aj, l)Q n (E, aj, I) , (n > 0, J = -(aj+2/+l)/4 G N) . 

(2-6) 

Thus, if Pj(E) vanishes then so do all P n >j(E) and the series (|2.4p terminates, 
automatically giving a normalisable solution to ()2.1|) . The J roots of Pj(E) are 
the J exactly-solvable energy levels for the model. For J = 1 and J = 2, the 
exactly-solvable eigenvalues are 

J = 1 : £ = 0, J = 2 : E± = ±2 v / 2\ / 3 + 2/. (2.7) 

For all real values of the parameter a and i > —1/2 the sextic eigenproblems 
(|2.1l 12. 2p have up to scaling exactly the same eigenvalues as a family of third-order 
eigenproblems [HI [12] . The relevant third-order differential equation is 



n 3 

where 



dx 3 x 3 Vx 2 cix x 3 



E <f>(x) (2.8) 



G = 2 - (g gi + g g 2 + 5152) , (2.9) 
L = -2- g gig 2 + (soffi + 5052 + , (2-10) 
50 + 51 + 52 = 3, (2.11) 



*In equations (|2.9[l and (|2.10p we have corrected an overall sign typo in the corresponding 
equations of [12] . 
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and the boundary conditions on the positive real axis are 



0U_>o = x 9l (l + 0(x 3 )) , (f>(x) = Oix' 1 e~ x2/2 ) as x -> oo (2.12) 

with go < gi < Q2- The asymptotic condition ensures the other two possible 
behaviours of the solution at infinity are ruled out. The isospectrality of T-L 2 , Hj, 
with boundary conditions (|2.2l 12. 12j) respectively was first discussed in [9] for 
I = and, with the help of results from [TJ], generalised to I 7^ in [12]. The 
result is that the eigenvalues E, E associated to H 2 , %3 satisfy 

E = E/k, k = 4/(3\/3), (2.13) 

whenever the parameters {ct,l} and {go, g±, g 2 } in the two models are related as 

a = 2(2 -go- g 2 ) , I = {2g 2 - 3 - 2g )/Q , (2.14) 

and 

5 o = (1- a- 6Z)/4, gi = (l + a/2) , g 2 = (7-a + 6Z)/4 . (2.15) 

The result was obtained by showing that the associated spectral determinants- 
functions constructed to vanish at the eigenvalues-are proportional. By analytical 
continuation from I to — 1— I, the isospectrality result also extends to / < —1/2 

IS 

If a is now tuned to the quasi-exactly solvable points a = — (4J + 2/ + 1) of H 2 
for positive integer J, then via the isospectrality (|2.13p the third-order problem 
W3 with 

g = l/2 + J-l, 51 = 1/2-2J-Z , g 2 = 2 + J+2l (2.16) 

has a hidden QES sector and J eigenvalues can be found exactly. Since the 
isospectrality proof relates the eigenvalues and makes no conclusions about the 
eigenf unctions, we cannot immediately state that H3 at the points (|2.16p is itself 
'fully' quasi-exactly solvable, if for full quasi-exact solvability one would insist on 
being able to find algebraically not only a subset of the eigenvalues but also the 
corresponding eigenfunctions. One might expect that H3 acts invariantly in a 
finite-dimensional subspace spanned by functions of the form x 91 exp(— x 2 /2) x n . 
If that is the case, the QES eigenfunctions will take a simple factorised form 
generalising (|2.4p . We shall show that this is not in general the case. 

In [13] and [14], the spectral link between % 2 and H3 was naturally extended 
to the adjoint operator %\ to (|2.8[) : 

X \x) = -E X ] (x) (2.17) 

where = G , V< = —L. The appropriate boundary conditions are specified by 
X f U^ = ^(l+O(x 3 )) , X f Uoo^0 , (gl = 2-g ) . (2.18) 



n\ xH 



x A + 

xy 



-dx? 



1 a 



1 



.T 2 dx, .T 3 
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It is important to note that the eigenfunctions <p n of H3 and the eigenfunctions 
Xn of Tl\ with eigenvalues E n and — E n respectively are substantially different 
functions, not related to each other by simple conjugation. Provided the roots of 
the indicial equation are ordered as go < gi < 52 then for g\ > — 1/2 the </» n 's are 
square integrable on M + , 

lim xl{x)<t>m{x) = O(x 2 " 90+91 ) ^ , lim xU^Mx) , (2.19) 

x— >0 x— >oo 

and the sets {xli} and {4> n } can always be normalized such that 

POO 

(n\m) = / Xn{x)(j) m {x) dx = 8 n , m . (2.20) 
J 

The properties described above are reminiscent of the well-studied properties of 
non-self-adjoint spectral problems and biorthogonal systems in quantum mechan- 
ics. These systems were introduced and studied in the early days of quantum 
mechanics, and were more recently revisited in the context of "PT-symmetric 
quantum mechanical models [15]. The reader is addressed to [16] for a recent 
review of this material. 



3 Bender-Dunne polynomials and projective triviality 

It is natural to ask whether a simple factorisation similar to (|2.4j) also characterises 
the exactly-solvable energy levels of and T-L\. To answer this question we first 
identify a simple necessary condition for the existence of a factorisable solution 
to equation (|2.ip of the form 

^(x) = x l+1 V 2 j-2(x,E)e- x4 / A , (3.1) 

where T , 2j(x,E) is a polynomial of order 2J in x. The large- £ behaviour of ij; 
should match the general WKB result 

V>(z) = 0(x~ 3/2 ~ a/2 e- x4/4 ) x -> 00 . (3.2) 

The ansatz (|3.ip agrees with (|3.2p for a = aj = —{21 + 1 + 47), a result which 
precisely matches the set (12. 3p . Let us now assume the existence of a solution to 
(1231) of the form 

( f>(x) = x^V 2 j-2(x,E)e- x ^ 2 . (3.3) 
The relevant WKB asymptotic for the third-order ODE (|2.8p is 

<j>(x) = 0{x- l e- x2/2 ) x -»• 00 . (3.4) 

Comparing with f|3.3|) we find 51 = 1 — 2 J, a very restricted set of values compared 
to (12. 3p . Starting from equation (12.17P leads to a similar conclusion. It is easy to 
check that when g\ = 1—2 J the ansatz (|3.3p produces a single eigenfunction with 
eigenvalue E = for all odd integers g\ satisfying go < gi < g2- However, away 
from these points we conclude that the eigenfunctions do not take the simple 
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factorised form (I3.3p . Although other types of wavefunction factorisation cannot 
be ruled out by this simple argument, we have checked that Tis is not one of the 
higher-order operators that appear in Turbiner's classification [3]. The surprise 
is that despite this, a subset of the eigenvalues of H3 with parameters (|2.16p can 
be found algebraically as zeros of certain polynomials. 

Since we may not have identified the QES sector of without the link to the 
isospectral problem H2 , it is natural to ask whether there are any alternative ways 
to detect the appearance of this hidden quasi-exact solvability in our third-order 
equations or in other models. 

To answer this question it is convenient to step back to [12] where, amongst 
other results, a series of full and partial isospectralities for T-L2(ot,l) = %2 with 
boundary conditions (12. 2|) were observed. Four of them are summarized by the 
following diagram 

^ 2 (-(4J+2/+l),/) — > n 2 (2J-2l-l,J+l) 



H 2 (2J+4l+2,-J-±) — > U 2 {2J+4l+2, J-\) 

Vertical arrows correspond to eigenproblems that have exactly the same eigenval- 
ues, while a horizontal arrow connects two problems that have the same eigenval- 
ues up to the elimination of all the QES levels present in the left hand models. 
The two problems on the bottom row correspond to the same Schrodinger equa- 
tion, and differ only in the boundary condition imposed at the origin. It follows 
from the diagram that the 'regular' eigenvalue problem for this equation, that 
with the behaviour x J+l / 2 at the origin, has exactly the same spectrum as the 
irregular problem with the x~ J+l l 2 behaviour at the origin, with the exception 
of the first J eigenvalues. 

For general I, the eigenproblem % 2 {2J+4l+2, — J— |) also does not appear in 
Turbiner's list of QES models [TJ. Nonetheless, it is isospectral to the QES sextic 
Schrodinger problem and J of its eigenvalues can be found exactly. In [12], after 
noticing an interesting symmetry in the recursion relation for the Bender-Dunne 
polynomials, it was (erroneously) stated that the appearance of QES eigenvalues 
in "%2(2J+4Z+2, — J— ^) corresponds also to a factorisation of the eigenfunctions 
in the form (|2.4p . The latter statement can be checked using the simple consis- 
tency criteria introduced above. Setting 

= x~ J+l ' 2 V 2K -2(x, E)e~ xA l A , (3.5) 

at large x the wavefunction behaves as x ~ ,J ^/ 2+2K e~ x4 / 4 while the WKB pre- 
diction is x ~^l 2 ~ J ~ 2l e~ x / 4 . Hence (13.51) is a suitable ansatz only for 

£ + 1/2 = -if , K = 1,2,3... . (3.6) 

Again, the constraint (|3.6p is much stronger than f|2.3|) and we should conclude 
that for general values of I the QES wavefunctions of %2(2J+4Z+2, — J— ^) do 
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not take a factorised form such as (|3.5p . We have not ruled out that the eigen- 
functions can be written in terms of other elementary functions, a point to which 
we return at the end of this section. However, the key point is that with the stan- 
dard techniques the quasi-exact solvability of %2(2J+4Z+2, — J— i) would not be 
evident and we would not know a-priori how to determine the QES eigenvalues. 

Returning to the question of how to detect quasi-exact solvability in such 
problems, we note that the eigenfunction ip(x,E,l) selected by the boundary 
conditions (|2.2p is one of two solutions to (|2.ip . characterised by their small- x 
behaviour 

if){x,E,l) = 0(x l+1 ) , ip(x,E,-l-l) = 0(x~ l ) x^O. (3.7) 

Provided the roots I and — l—l of the indicial equation do not differ by an integer, 
then the ^-functions (I3.7P are linearly independent. Moreover the solutions are 
automatically projectively trivial around the origin, by which we mean that for 
arbitrary E the monodromy of ip(x, E, I) around x = is such that 

i>(e i2w x, E, I) = e i2 < l+1 ^(x, E, I) . (3.8) 

The monodromy of if)(x,E, —l—l) follows via analytical continuation I — > —l—l. 

When I = — J — \ the roots of the indicial equation differ by 2 J and so for J 
integer there is a 'resonance' between the two solutions ip(x, J) = ip(x, E, — J— |) 
and tJ)(x,—J) = ip(x, E, J— i). This pair is no longer a basis of solutions to %2 
and a linearly independent solution i/)(x, E) to ip(x, E, J) must be constructed. In 
contrast to the regular solution ip{x,—J), the final solution will in general posses 
an additional logarithmic branch point at x = and thus the projective triviality 
property will therefore be lost: 

^(e i2w x, E) + e i2 ^- J+1/2) i>(x, E) . (3.9) 

We now show that the eigenfunctions corresponding to the QES eigenvalues of 
H 2 (2J+Al+2, - J-\) do not acquire logarithmic terms and so remain projectively- 
trivial. We suggest that this may be a means of detecting the hidden quasi-exact 
solvability of such models. To illustrate this we apply the Bender-Dunne method 
to n 2 (2J+4l+2, -J-\) by setting 

i>(x) = e- 4 / 4 x- J+1 ' 2 f; {-\) n 9AErIA x 2n (3 1Q) 

n=0 

In comparison with (|2.4p . the gamma function has been dropped from the de- 
nominator of ()3.10p to ensure the coefficients in the series remain finite for all n 
and (I3.10P is a -by construction- projectively-trivial solution ip(x). Consequently 
the polynomials Q n (E) now satisfy 

(n - J)Q n (E) = EQ n -x(E) + 16(n - l)(n + I - l/2)Q n _ 2 (£) . (3.11) 

Setting Qo(E) = 1, the recursion relation defines Q n (E) in terms of E and J for 
n < J just as before. The first difference occurs when n = J since the LHS of 
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(|3.1ip vanishes. The RHS is a J th -order polynomial in E which must be zero. 
Its roots therefore determine the QES eigenvalues E n , n = 1, ... J. However, the 
recursion relation (|3.1ip has left Qj{E) unspecified and provided I is such that 
n+l— 1/2 7^ the remaining coefficients do not factorise. In general, all Q n >j{E) 
will be a function of the unknown coefficient Qj(E). For example 

Qj+i(E) = EQj(E) + 16J(J + I + l/2)Qj_i(£) . (3.12) 

Since Qj-%(E n ) ^ 0, from (13.12P we see that the choice Qj(E n ) = does not 
lead to a truncation of the series. More precisely, the resulting solution (|3.10[) will 
always have the desired monodromy properties but only for a precise value of the 
constant Qj(E n ) will it be asymptotically subdominant and satisfy the boundary 
condition at infinity. 

Returning to the question of determining the exact eigenfunctions of the fam- 
ily %2(2</+4Z+2, — J— 1/2), we note that for all odd integers J one of the QES 
eigenvalues is E = 0. In this case the eigenproblem can be solved in terms of a 
Whittaker function 

IK*) = 2J/4r / !+ 3 | + " ) W j ! ufe) , (3.13) 
which exactly matches (|3.10p when E = for the choice 



(_i)j/2-i/2 2 5J / 2 y^r(f+i ) r(f +f +|) 
r«)r(i+i) 



Noticing that for small J the solution (|3.13p when I is an integer can be written 
in terms of Bessel functions multiplied by polynomials in x, we were motivated 
to try an ansatz of the form 



= X ^-J (a n (E) x 2n K L (f ) + b n (E) x 2n kJ^U (3.15) 



4 V ^ / 4 
n=0 



where K n {x) are modified Bessel functions of the second kind and a n (E),b n (E) 
depend on E, I. Acting with H 2 (2 J+4/+2, - J- 1/2) for integer I on (^T5|) and 
comparing powers of x, we found for J = 1 that E must be zero, the series on the 
RHS of (|3.15p truncated at n = 1 + 1 and (|3.15p reproduces (|3. 13|) as anticipated. 
Setting J = 3, we solved for the coefficients {oq, «i, . . . , a^ + i, bo, b\, . . . , 63+/} and 
found one solution with E = and two further solutions with eigenvalues E = 
±8v2 + I, exactly reproducing the three solutions of Ps(E) = 0. Repeating this 
process for odd J, we found that (|3.15p generates J wavefunctions and constrains 
the QES eigenvalues to be solutions of Pj(E)=Q. The series on the RHS of (I3.15P 
truncated at n = J + I. 

Given the ansatz (]3. 15|) unexpectedly generated not just the zero eigenvalues 
but all of the QES eigenvalues for I integer when J is an odd integer, we then 
checked if it also works for even integers J. It turns out the ansatz (13.15P indeed 
gave two solutions satisfying the boundary conditions (|2.2p when J = 2 provided 
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E = ±2\/2V3 + 21, exactly matching ([221). We find that the ansatz (j3T5l) worked 
for all integers J, I £ N + with the series on the RHS of (|3.15p truncating at 
n = 2(J + I). However, at present we are not able to generalise this ansatz to 
non- integer values of I. 

Differential operators that act invariantly on a subspace spanned by poly- 
nomials multiplied by special functions of either hypergeometric, Airy or Bessel 
type have been constructed in [17] . With the variable change x = \^2w 1 ^ 4: and 
the gauge transformation w~ 9 / 8 %2 w 5 / 8 , we find 7i2(2J+l+2, — J— 1/2) becomes 
proportional up to an additive constant to the differential operator presented 
in [T7] for the cases when E = and I is an integer. To reproduce the solutions 
(|3.13p . the invariant subspace given in [T7] for J§ must be extended to include 
Bessel functions multiplied by rational powers of x. We leave further details of 
these exact wavefunctions and the investigation of the cases when I is not an 
integer to future work. 

4 Projective triviality and third-order QES models 

Continuing the discussion of the last section, we now show that the hidden QES 
sectors of H3 and can be detected by using the projective-triviality test dis- 
cussed above on x^ = x[ )> one °f the three linearly independent solutions {xIa} 
to (|2.17p . These solutions are characterised by their small- X behaviour 

X \ i) (x,E) = 0(x°t) , x^O (4.1) 

where g\ = 2—§i(i = 0, 1, 2) are the roots of the indicial equation. These solutions 
are projectively trivial around the origin 

x\ l) (e^x,E) = e^xl l) (x,E) (4.2) 

for (gi — gj) ^ Z with i / j. The QES eigenvalues of ^3, H\ appear when 

go = 1/2 + J —I , gi = 1/2 -2J-1 , g 2 = 2 + J + 21 , (4.3) 

and 

gl = 3/2 - J + I , g\ = 3/2 + 2J + I , g\ = — J — 21 . (4.4) 
Hence we have 

gl=g\-3J , J = 1,2,3... . (4.5) 

Furthermore, for / > —1/2 and J > the ordering is g\ < g\ < g\, and we are 
in the presence of a resonance phenomena for the solutions {xj}- Again, this 
circumstance usually leads to the loss of the projective triviality property (14. 2p 
due to the appearance of logarithmic contributions to the wavefunction. The 
analysis of $3] suggests that imposing projective triviality on the eigenfunctions 
may be a way to identify the set of exactly known eigenvalues. 
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Instead of using a Bender-Dunne like ansatz for the wavefunction we will 
construct x perturbatively using Cheng's method [18J. The solution x to (|2.17[) 
is also a solution to the equation 

x \ x ) = x 2 ~9o + L [( x 3 _ E) X ] {x)\ , (4.6) 

where 

L ^ = n ( TTJ S ' (A; = 0, 1, 2) . (4.7) 
Uk(P+ 1 + 9k) 

The function %t can be considered as the n — > oo limit of a function x^ n ^ con " 
structed from x^ "* = x 2_So using the following recursion relation 

x t(n) (x) = x t(0) (x) + L[(x 3 _ £) x t(»-l) (a .)] . ( 4. 8) 

After a single iteration, we have 

y t(i) = x 2- ff0 ^! + x l "I f 49 ) 

V Ylk(Z-90+9k) 11^-90 + 9k) ) ' K '> 

Using (j4.9|) we can study the QES problem at go — gi = 3 J when J = 1. As 
(<7o — 5i ) —^3 the second term in the parenthesis on the RHS of (|4.9p is in general 
divergent. However if the limit (go — g\) — > 3 is taken simultaneously with E — > 
so that E/(3 — go + gi) — > C with C finite, the final result is again finite and, 
up to this order in the perturbative expansion, x remains projectively trivial. 
The result E = P\{E) = corresponds precisely to the only exactly-solvable 
energy level at J = 1. Further, the result of a second iteration is 



x t(2) = x 2- 90[1 



Ex 3 



x 
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]lfc( 3 -9o + 9k) n fc (6 -go +gk) 

Ex 6 E 2 x % 
+ 



llfc( 3 -9o + 5fc)(9 - go + gk) ]lfc( 3 -9o+ 5fc)(6 - go + gk) 

Ex 9 x 12 , 



llfc( 6 - 50 + 9fc)(9 - 5o + £fc) llfc( 6 - 90 + 5fc)(12 - 5o + 0fc) 

At J = 1, the potential divergences again disappear in the limit Ej (3 — go + g\) — > 
C, and it is easy to check that simultaneously all the subsequent x^ n ^ s remain 
finite. It is also possible to check that the solution corresponding to C = does 
not lead to the desired subdominant solution (I4.6p . The exact solution for C = 
is 

X t (x)|c=o = x 2 -^ F 2 [i,2-f,^] , (J =1,^ = 0) (4.11) 

which indeed grows exponentially as x~ 1 e x2 ^ 2 . In order to find the proper wave- 
function we use the asymptotics 

z (2-9o)/6 oF2[ l 2-%,z]~ r(2 ~jV l/6 e 3^ 3 z ; (4 . 12) 

2v/ 3tt 
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and introduce the linearly independent solution 

x\ = z {5 - go)/ %F 2 [ll-f,z} (4.13) 
which behaves asymptotically as 

p/5 _ 9o\ 

X \ ~ ik^i z -i/6 e ^ ^oo. (4.14) 
4v37r 



Thus, the asymptotically vanishing solution is 

/ p/o _ 90 ^ 

t(V) _ x 2-go n i? 9 [I 2-20 «1] 2j_ 3 p r3 5 _ 90 

A W X I 0^2 L 2 ' ^ 2 ' 21fiJ o /^WK 90^ 2 ^2' 2 2 > 21f 

T 



2 ' 216 J 3 y /g|-(5 _ _ " - L 2 " 2 2-21.0-1 



(4.15) 

The solution flng) corresponds to the choice C = 3 v / 6r(2 - f )/r(| - f ). 

The case J = 2 can be treated in a similar fashion: as (go — g\) —> 6 the 
x e coefficient of (I4,10p diverges. This singular behaviour can be avoided in the 
double limit 

_ _ lim (l+o7o — 4^ —^)/(G-9o + 9i)^C (4.16) 

e-*e ± , (90-91)^6 V 3 ( 3 - 9o +ffi)(3 - 90 + 92) J 

provided 



£± = ±3\/3a/-2 - 51 = ±3y - V3 + 2/ . (4.17) 

The result (|4.17p matches the exact energy levels (|2.7p provided E/E = k = 
4/(3\/3) (cf. (|2.1ip ). Further, the numerator of (|4.16p is simply related to the 
Bender-Dunne polynomial P2(E) (see below) and it is possible to argue that there 
always exists a value of C such that the wavefunction decays exponentially at large 
x. This proves, therefore, that the roots of the Bender-Dunne polynomials P%(E) 
are indeed part of the spectrum of the dual pair (12. 8( I2.17P . 

More generally, for any {gi} the full Cheng solution can be written in the form 

Vn^O rifc=o( 3n -50+5fc)/ 

where P n (E) are degree n polynomials in E that satisfy 

2 

P n (£) = EP n ^(E) + JJ(3(n - 1) - g + g k )P n -2{E) (4.19) 

fc=0 

with Pq = 1,Pi = E. Restricting {gi} to the QES points f|4.3[> . the recursion re- 
lation (|4.19p matches the Bender-Dunne recursion relation (|2.5p with n n P n (E) = 

Pn(E). 
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5 Other models and conclusions 



In this paper we have reported some progress toward a more complete understand- 
ing of the spectral equivalence between the sextic anharmonic oscillator (12. If) and 
the dual pair of third-order ODEs (|2,8p and (|2,17p , and the nature of the quasi- 
exact solvability of the third-order problems that is thereby induced. Although 
more work will be needed to complete the picture, and more generally to un- 
derstand the emergence of standard and hidden quasi-exact solvability in dual 
pairs {'H,'H^} of higher-order differential operators, we think that the concept of 
projective triviality should be a useful tool in detecting QES sectors. In addition 
to the examples discussed in £}3] and $H we have discovered that hidden quasi- 
exact solvability is a property shared by many other models. For example, we 
have applied the same analysis to the n th order differential equations introduced 
in [19] 

[(-l) n+1 D(g n ^ - (n-1)) D(g n ^ 2 - (ra-2)) . . . D(g ) + x nM ] ^(x) = E^x) 

(5.20) 

where 

°(^(|-f) • E» = =^ <-> 

and the boundary conditions are 

VUo = O(x^) , iP = 0(x^- n ^ 2 e - M+1 /(A/+D) as x ^ oo (5.22) 

with go < gi < • • • < <7n-i- These directly generalise the problems 7^2(0,— go) 
and %3. By imposing projective triviality on the wavefunctions of the adjoint 
problems at the resonant points go — gi = nJ for M, J G N, we found the Cheng 
solutions for the adjoint problems are 

*•<*) - ^- ( £ u J-^ iE)x Z ) <»■»> 

v^o n J= i n fc =o w -9o+ 9k) j 

where the corresponding Bender-Dunne polynomials satisfy 

M n-1 

P m (E) = EP m _i(E) - (-1) M 11 JJ( n (j+m-M-l) - g Q + g k )P m -i- M (E) . 

j=l k=0 

(5.24) 

The QES eigenvalues are the J roots of Pj(E) = and, in general, the associated 
wavefunctions do not have an elementary form. As noted in for the second- 
order models, the non-QES part of the spectrum is precisely the spectrum of the 
same differential equation subject to a boundary condition that imposes regular 
behaviour of the wavefunction at the origin. 

Finally, we should reiterate that for all the QES models encountered in this 
paper when go — gi = nJ, the ordering go < 9i < 92 < ■ ■ ■ of the solutions of the 
indicial equations is not fulfilled. Hence, the norm (12,201) of the exactly-solvable 
states \J (n\n) is divergent and the corresponding radial eigenvalue problem is 
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always 'irregular'. Higher-order differential equations of the form (|5.20p have 
recently been studied in the context of VT symmetric quantum mechanics and its 
generalisations for even n in [20} I2T] . Additional motivation for the further study 
of higher-order eigenproblems of the type considered in this paper comes from 
their relevance to particular integrable quantum field theories, via the ODE/IM 
correspondence [221 1231 [T9l H31 123] . 
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